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Abstract. The turbulent transport of angular momen- 
tum plays an important role in many astrophysical ob- 
jects, but its modelization is still far from satisfactory. We 
discuss here what can be learned from laboratory experi- 
ments. We analyze the results obtained by Wendt (1933) 
and Taylor (1936) on the classical Couette- Taylor flow, 
in the case where angular momentum increases with dis- 
tance from the rotation axis, which is the most interesting 
for astrophysical applications. We show that when the gap 
between the coaxial cylinders is wide enough, the criterion 
for the onset of the finite amplitude instability can be ex- 
pressed in terms of a gradient Reynolds number. Based 
on Wendt's results, we argue that turbulence may be sus- 
tained by differential rotation when the angular velocity 
decreases outward, as in keplerian flows. From the rota- 
tion profiles and the torque measurements we deduce a 
prescription for the turbulent viscosity which is indepen- 
dent of gap width; with some caution it may be applied 
to stellar interiors and to accretion disks. 
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1. Introduction 

Differential rotation is observed in various astrophysical 
objects, from planets to galaxies, and one suspects that 
it gives rise to turbulence, since shear flows are liable to 
hydrodynamical and MHD instabilities. When these in- 
stabilities are of the linear type, they are relatively easy 
to study by perturbing slightly the equilibrium state. But 
some of them occur only at finite amplitude, in which case 
the answer must be sought in computer simulations or lab- 
oratory experiments, with their inherent limitations. 

There has been some debate recently on whether a ke- 
plerian disk, which is linearily stable (Rayleigh 1916), may 
be unstable to finite amplitude perturbations. It may look 
as if this question presents little interest, since it has been 
proved that a very weak magnetic field suffices to render 
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such a disk linearly unstable (Chandrasekhar 1960; Balbus 
and Hawley 1991). However the properties of angular mo- 
mentum transport depend sensitively on which instability 
dominates in the considered regime, and a finite ampli- 
tude instability can overpower the linear instability which 
is the first to occur, as the relevant control parameter in- 
creases. One example is the Couette- Taylor flow, with the 
outer cylinder at rest and the inner cylinder rotating with 
angular velocity fl. When £1 is increased, the transport 
of angular momentum first scales as f2 3 / 2 , but thereafter 
it varies as ft 2 , once the flow has become fully turbulent 
(Taylor 1936), as if the initial linear instability were su- 
perseded by a stronger shear instability (see also Lathrop 
et al. 1992). By extrapolating the O 3 / 2 law to high Q one 
would clearly underestimate the transport. 

In the present article, we take as working assumption 
that any differentially rotating flow experiences, at high 
Reynolds number, the turbulent regime observed in the 
Couette- Taylor (CT) experiment, and that this turbulence 
will then transport angular momentum in the same way as 
in that experiment. The CT flow has been chosen as refer- 
ence because it is the simplest flow to realize in the labo- 
ratory, with both shear and rotation that can be varied in- 
dependently. We examine whether the experimental data 
suggest a prescription for the angular momentum trans- 
port, which may be used to model astrophysical objects. 
A similar approach has been taken by Zeldovich (1981), 
but our conclusions will differ from his (see Appendix). 



2. The Couette- Taylor experiment 

The CT experimental apparatus consists of two coaxial 
rotating cylinders of radius R\ and R2 separated by a gap 
Ai? = i?2 — Ri, which is filled with a fluid of viscosity v. 
The cylinders can rotate with different angular velocities 
Vt\ and 0,2] their height is in general much larger than 
their radius, to minimize the effect of the boundaries. The 
Reynolds number is usually defined in terms of the dif- 



2 



Denis Richard and Jean-Paul Zahn: Turbulence in differentially rotating flows 



5.00 



4.50 - 



/ 

L 

/<o- 

/ »" 

•/ 



« 4.00 

a 

< 



3.50 



• / 

/ 

/ 

/ 



/ Re = 2000 

(plane shear flow ) 



3.00 



-2.00 



-1.50 -1.00 
log(AR/R) 



-0.50 



0.00 



Fig. 1. Critical Reynolds number vs. aspect ratio; 
filled circles from Taylor (1936), open circles from 
Wendt (1933). Dotted line: critical Reynolds number for 
plane shear flow instability; dashed line: critical gradient 
Reynolds Re* for circular shear flow in the limit of large 
gap (see text). 



ferential rotation AS1 = | O2 — fii | and by taking the gap 
width as characteristic length: 



Re 



AilRAR 



(1) 



where R is the mean radius R = (i?i + i?2)/2. 

2.1. Critical parameters and transition to turbulence 

When the inner cylinder is rotating and the outer one 
is at rest, angular momentum decreases outward and the 
flow is linearly unstable for Reynolds numbers higher than 
Re c = 41.2 ^/R/AR (for narrow gap, cf. Taylor 1923). 
This well known instability takes first the shape of steady 
toroidal, axisymmetric cells (the Taylor vortices); it is very 
efficient in transporting angular momentum, whose gradi- 
ent is strongly reduced. At higher Re a wavy pattern ap- 
pears and after a series of bifurcations the flow becomes 
fully turbulent. This case has been studied by many exper- 
imental teams, and it is extremely well documented (see 
Andcreck et al. 1986). It has also been modeled success- 



fully in three-dimensional numerical simulations (Marcus 
1984; Coughlin & Marcus 1996). 

In the opposite case, when the outer cylinder is rotat- 
ing and the inner one is at rest, the angular momentum 
increases outward and the flow is linearly stable. The only 
theoretical prediction concerning the non-linear behavior 
is that by Serrin (1959), later refined by Joseph and Mun- 
son (1970), who established that the flow is stable against 
finite amplitude perturbations below Re = 2n 2 (in the 
narrow gap limit). 

To this date, no numerical simulation has been able 
to demonstrate the finite amplitude instability. But this 
instability does occur in the laboratory, and it has been 
described already by Couette (1890). It was studied in de- 
tail by Wendt (1933) and Taylor (1936), who showed that 
for Reynolds numbers exceeding Re c ~ 2 10 3 , the flow be- 
comes unstable and immediately displays turbulent mo- 
tions. The critical Reynolds number depends on whether 
the angular speed is increased or decreased in the exper- 
iment, a typical property of finite amplitude instabilites. 
Moreover, it is sensitive to gap width, as demonstrated 
by Taylor. Figure 1 displays results from Wendt and Tay- 
lor: Re c is roughly constant below AR/R = 1/20, but 
above it increases as (AR/R) 2 , as was already noticed 
by Zeldovich (1981), a behavior for which an explanation 
was proposed by Dubrullc (1993). In the latter regime one 
can define another critical Reynolds number Re* involv- 
ing, instead of gap width, the gradient of angular velocity; 



since 



Re, 



R? Aft ( AR\ ' 
~V~AR \~R~) 



the instability condition becomes 



Re* = — 



R 3 An 



AR 



> Re* ~ 6 10 5 



(2) 



(3) 



We see that two conditions must be satisfied for the finite 
amplitude instability to occur: the first Re > Re c is the 
classical criterion of shear instability, valid also for plane 
parallel flows, whereas the second Re* > Re*, involving 
what we shall call the gradient Reynolds number, is gen- 
uine to differential rotation. In addition, to trigger the 
instability the strength of the perturbation must exceed a 
certain threshold, which presumably also depends on Re 
or Re*. 

2.2. Transport of angular momentum 

In the turbulent regime, the torque measured by Tay- 
lor scales approximately as G a (^2)™ for a given gap 
width, where the exponent n tends to 2 for large Q2 ■ The 
measurements made by Wendt confirm that scaling with 
n w 2. It suggests that the transport of angular momen- 
tum may be considered as a diffusive process, and that the 
mean turbulent viscosity V t increases linearly with Q 2 > or 
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Fig. 2. Value of parameter a in the classical viscosity 
prescription v t = aRAil AR, derived from Taylor (1936) 
in filled circles and from Wendt (1933) in open circles. 

AO. It is then natural to examine whether this viscosity 
may be expressed as 



v t = aR An AR 



(4) 



with a being a constant of order unity, since the largest 
turbulent eddies would have a size « AR and a periph- 
eral velocity « R AO/2. The parameter a is easily derived 
from the torque measurements, and the surprising result 
is that it decreases with gap width (fig. 2). For the small- 
est gaps, a scales as the inverse of AR/R, but the slope 
steepens farther as if the scaling would tend asymptoti- 
cally to 



a oc 



(¥)' 



f ( AR 



i. 



(Comparing Wendt's experimental data of his two largest 
gaps, one deduces an exponent 6\na/Sln(AR/ R2) = 
-1.83.) 

We may therefore conclude that, in the limit of large 
gap, the mean turbulent viscosity actually scales as 



aAR\ 3 AO AO 
-R-J R AR =f3R AR' 



(6) 



with (5* w 4 10" 



This strongly suggests that the local value of the turbu- 
lent viscosity is then independent of gap width, and that 
it is determined only by the local shear: 



v t = f3r 3 



(19. 



dr 



(7) 



r being the radial coordinate. 

In principle, one should be able to verify this prescrip- 
tion for v t by examining the rotation profiles measured by 
Taylor and Wendt. According to (^), the conservation of 
angular momentum requires that its flux, given by 



T 



f3r 3 



dr 



,dO 

dr ' 



(8) 



varies as 1/r between the cylinders. Therefore r 3 (dfl/dr) 
should be constant in the turbulent part of the profile (as 
it is in the laminar flow). 

But this constancy can be expected only if the trans- 
port of angular momentum is achieved by the viscous and 
turbulent stresses alone. That is not the case in Taylor's 
experiment: as acknowledged by him, an Ekman circula- 
tion is induced by the ends of his apparatus, although he 
tries to minimize the boundary effects by chosing a large 
aspect ratio (height/radius). Moreover his rotation pro- 
files arc deduced from pressure measurements made with 
a Pitot tube located at half height of the cylinders, where 
the radial return flow has its maximum intensity. Conse- 
quently, the torque inferred from these profiles is about 
half of that measured directly at the inner cylinder. 

The aspect ratio is less favorable in Wendt's experi- 
ment, but there the top boundary is a free surface, and 
most of his measurements have been made with the bot- 
tom boundary split in two annuli, attached respectively 
to the inner and the outer cylinders, which reduces drasti- 
cally the circulation and renders his results more reliable. 
We examined his rotation profiles obtained for the largest 
gap width (AR/R = 0.38) and with 4 different speeds of 
the inner cylinder (0 < Oi < O2/2); we found that in the 
bulk of the fluid, these profiles are compatible with the 
constancy of the gradient Reynolds number, as predicted 
by (||). However we cannot rule out a mild variation of (3 
within the profile. Unfortunately, Wendt gives the torque 
only for the case where the inner cylinder is at rest; we 
draw from it the following value of the coefficient f3: 



( 5 ) = 1.5 ±0.5 IO- 



CS) 



The estimated uncertainty reflects that of the velocity 
measurements: /3 results from second derivation of the 
shape of the fluid surface. 

3. Are keplerian flows unstable? 

More precisely, the question we address is whether the 
Couette flow is unstable when the angular velocity de- 
creases outwards and the angular momentum increases 
outwards, as in keplerian flow. 
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Fig. 3. Angular velocity and momentum profiles in the 
case of decreasing angular velocity and initially constant 
angular momentum - dotted line (experimental data from 
Wendt 1933). 
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There is no definite answer yet, because this regime has 
not been explored at high enough Reynolds number. But 
some information can be gleaned from Wendt 's study. He 
reports the results of experiments where the two cylinders 
rotate such that ^2-^2 = ^ii£?- At low Reynolds number, 
this setup enforces a laminar flow of constant angular mo- 
mentum (neutral flow) , but at high Reynolds number this 
flow becomes turbulent for two of the three gaps used by 
Wendt. 

The angular velocity and angular momentum profiles 
for one of these turbulent flows are reproduced in fig. 3. 
(According to Wendt's data, f^-Ri actually exceeds fii-Rf 
by a half a percent.) The profiles clearly demonstrate 
the flow instability, with angular momentum being trans- 
ported down the angular velocity gradient, which becomes 
somewhat hater far enough from the boundaries, whereas 
the angular momentum profile steepens substantially. 

In the turbulent bulk of the flow q = —dlnd/dhir s» 
1.4, compared to the initial q — 2. We recall that q = 1.5 
in keplerian flow, and that in the numerical simulations 
performed by Balbus et al. (1996, 1998), the instability is 
lost already at about q — 1.95. A crude estimate of the 
parameter (3 in the viscosity formulation (Q) indicates that 
the size of the turbulent eddies is much smaller than the 
gap width: I = v7?i? w Ai?/100. 

The corresponding values of Re* for the three experi- 
ments are reported in fig. 4 together with the critical line of 
fig.l. They are located respectively above and below this 
line for the unstable and stable flows. The data are too 
scarce to locate precisely the critical line of these "neu- 
tral" flows, but we can conclude that the critical gradient 
Reynolds number Re* then lies between 2 10 5 and 6 10 5 . 



Fig. 4. Reynolds number vs. aspect ratio for the three 
"neutral" experiments from Wendt: the two filled circles 
correspond to the unstable flows, the open circle to the sta- 
ble flow. The stability curve obtained for the inner cylinder 
at rest (fig.l) is displayed for comparison. 

4. Discussion 

The firmest result of our analysis of the Couette- Taylor 
experiment is that the criterion for finite amplitude in- 
stability may be expressed in terms of the gradient 
Reynolds number Re* — R 3 (Afl/ AR)/v, where the crit- 
ical Reynolds number Re* c < 6 10 5 is independent of 
the width of the gap between cylinders, for wide enough 
gap. The turbulent transport of angular momentum then 
seems also to be independent of gap width; it proceeds 
always down the angular velocity gradient, as confirmed 
by the behavior of the initially "neutral" flows examined 
by Wendt. 

Though the experimental evidence is somewhat less 
compelling, we have established empirically an expression 
which links the turbulent viscosity to the local shear. The 
value of i?e*, and that of f3 in eq.(0), have been derived 
from Wendt's experiment with the inner cylinder at rest, 
and it is not obvious that these parameters would be the 
same for different ratios of cylinder speeds. Also, the linear 
scaling v t jv — (3Re* may be valid only for those moderate 
gradient Reynolds numbers which could be reached in the 
laboratory 

Nevertheless, it is tempting to apply this expression (j7|) 
to accretion disks, as an alternate for the commonly used 
prescription v t = ac s H, where H is the scale height of the 



Denis Richard and Jean-Paul Zahn: Turbulence in differentially rotating flows 



5 



disk and c s the local sound speed (Shakura & Sunyaev 
1973). Some caution is required because this viscosity has 
been derived from experiments performed with an incom- 
pressible fluid. Moreover (Q) implies that the eddies which 
dominate in the transport of angular momentum have a 
size of order I » \f]5 r, independent of the strength of the 
local shear, and that their velocity is of order £rd£l/dr. 
When applying this prescription to a compressible flow, 
one has to make sure that this velocity is smaller than the 
sound speed and, in the case of an accretion disk, that the 
size £ of the eddies, which are three-dimensional, does not 
largely exceed the scale height H . The behavior of "neu- 
tral" flows demonstrates that the shear instability always 
transports angular momentum down the angular velocity 
gradient, which means outward for accretion disks. 

Note that in a keplerian disk our expression is equiva- 
lent to 



with 



(10) 



Such a prescription has been suggested originally by 
Lynden-Bell and Pringle (1974), and recently it was used 
again by Duschl et al. (1998). As a test, it is being applied 
to the modelling of accretion discs in active galactic nuclei 
(Hure & Richard 1999). 

The reader may wonder why we have only used exper- 
imental results dating from the thirties, namely those of 
Wendt (1933) and Taylor (1936). The reason is that no 
one, since them, has studied in such extent the regime of 
outward increasing angular momentum.^ We suspect that 
it is because the flow becomes then turbulent at once, 
without undergoing a series of bifurcations associated with 
enticing patterns. But we hope that experimentalists will 
turn again to this classical problem, which is of such great 
interest for geophysical and astrophysical fluid dynamics, 
and that they will explore the rotation regimes for which 
the data are so incomplete. In the meanwhile, the quest 
will continue to detect the finite amplitude instability in 
computer simulations. 
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Appendix 

On Zeldovich 's analysis of Taylor's results. 

Zeldovich (1981) had a similar goal in mind when he 
analyzed the results of Taylor's experiment. But he started 
from the idea that the turbulent flow was governed by 
the epicyclic frequency Nq and the turnover frequency u, 
where 



Nh = ±-Mr^f 

r° dr 



and 



dfl\ 



since they measure respectively the stability of the 
flow and the strength of the shear. He defines a non- 
dimensional parameter Ty = N^/uj 2 , akin to the Richard- 
son number used in stratified shear flow, and he seeks the 
confirmation of his intuition that Ty be constant in Tay- 
lor's turbulent rotation profiles, in which case the angular 
velocity would obey a power law 

n cx r q . 

His best fit yields q — 5.5. However these profiles are 
contaminated by the Ekman circulation mentioned above, 
and they differ marquedly from those obtained by Wendt. 

A consequence of this constant Ty would be that the 
parameter a in (Bj) would vary as 



2q+4 



IS 



a property which is not substantiated by the combined 
results of Taylor and Wendt, as can be seen in Fig. 2. 



